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[Document Name] SPECIFICATION 
[Title of the Invention] Method of ML^EM image reconstruction 

and medical imaging apparatus 

[Claims] 

1. A method of ML^EM image reconstruction feruse in connection with a diagnostic 
imag ing apparatus that generates projection data comprising: 

(a) collecting projection data including measured emission projection dale; 

(b) assuming an initial emission map; 

(c) fjffiwwing an initial attenuation map* 

(d) iteratively updating the emission map and the attenuation map; 

(e) with each iteration, recalculating the emission map by taking a previous 
emission map and adjusting it based upon: © the measured emission projection data; (H) 
areprojection of the previous emission map which is carried out with a multi-dimensional 
projection mo del; and Q£) a reprojection of the attenuation map; and 

(f) with each iteration, recalculating the attenuation map by taking a previous 
attenuation map and sdjusmig h based upc^ ©them 

<H). the reprojection of the previous emission map which is carried out with the multi- 
dimensional projection model 

2. A method as claimed in claim 1, wherein the step of collecting projection data 
further includes collecting measured transmission projection data. 

3 . A method ss dsimedmclsim 2, wherein the adjustment to me pilous attenuation 
map in the r ecalculating step is also made based upon the measured transmission projection 
data. 
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4. A method as claimed m daim2 or claim 3 r wherein the step of ^n^^ j p^j^tf^ 
data further includes collecting a blank get of transmission projection data which la 
collected without a subject present in the diagnostic imaging apparatus, and the adfrstmmt 
to the previous attenuation nap hi the r ec a lculat ing itep it also made baaed upon the blank 
set of transmission projection data. 

5. A method as darned inany one of claims 2 to 4, wherem the measured transmission 
projection data Is truncated. 

6. Medical imaging apparatus for reconstructing image representations of a subject 
being examined therewith comprising: an emission memory (112) which stores measured 
emission projection data from the subject; an image memory (140) which stores miss ion 
maps; an attenuation map memory (130) which stores attenuation maps; a first p rojec to r 
(ISO) which generates emission map projections via fbrward projection of the emission 
maps from the image memory (140), said projector (150) using a first muMHfimenatooal 
projection model; a second projector (152) which generates attenuation map projections via 
forward projection of the attenuation maps from the attenuation map memory (130), said 
projector (130) using a second multidimensional projection model; a first data processor 
(160) which samples the emission memory (112), and the first projector (1S0X and in 
accordance therewith updates the attenuation maps stored in the attenuation map memory 
(130); a second data processor (170) which samples the emission memory (112), the first 
projector (150), and the second projector (152), and in accordance therewith updates the 
emission maps stored in the image memory (140); and a human-viewable display (910) 
which renders reconstructed images representations of the subject from the emission maps. 
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7. Medical imaging apparatus ts claimed in dam 6, further comprising: a trsnrrmwm 
memory (114) which stores measured transmission projection data from the subject, 
wherein the medical diagnostic imaging apparatus (10) collects said transmission projection 
data. 

S. Medical imaging apparatus as claimed in claim 7, wherein the first data processor 
(160) also samples the transmission memory (114) when it updates the attenuation maps 
stored in the attenuation map memory (130), such that the update is also made in 
accordance therewith. 

9. Medical imaging apparatus as claimed in claim 7 or dairn 8, fiirther comprising: a 
reference memory (120) winch stores transmission projection data from a blank run 
performed without the subject present in the medical diagnostic imaging apparatus (10X 
wherein the first data processor (160) also samples the reference memory (120) when it 
updates the attenuation maps stored in the attenuation map memory (130), such that the 
update is also made in accordance therewith. 

10, Medical imaging apparatus as claimed in any one of claims 7 to 9, wherein the 
measured transmission projection data from the subjected stored in the tr ansmissi on 
memory (114) is truncated. 
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[Detailed Description of the Invention] 
[0001] 

[ Technical Field to which the Invention pertains] 
The present invention relates to imaging using image reconstruction of e mis s io n dm. It 
finds particular application in conjunction with single photon nmiintoP computed 
tomography (SPECT) and nuclear cameras, and will be described with particular reference 
thereto. However, it is to be appreciated that the present invention is also amenable to other 
like applications including, but not limited to, positron emission tomography (PET) and 

ultrasound tomography. 

[0002] 
[Prior Art] 

Attenuation correction is row an integral part of nuclear med i ci ne image reconstruction. 
Several techniques have been introduced over die yean, ranging from simple analytical 
filters to elaborate transmission devices, to obtain an object-dependent a tte n u a tion map. 
Often, fa analy tical correction, assumptions are made that are too restrictive to be realistic 
Object-dependent attenuation collection based on emission data only is an attractive and 
theoretically sound alternative, but not ready for large scale utilization. Finally, techniques 
baaed on an atte nuation map constructed from measured transmission projection data, while 
a popular approach, suffer fiom very poor statistical quality. Moreover, current techniques 
for reconstructing the emission data are sensitive to this noise. Ad dit i onally , the collected 
transmission projection data is often truncated. 
[0003] 

The approaches presented herein are a departure from previously available techniques In 
order to appreciate the difference, it is helpful to examine previous techniques. 
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h is usually recognized that iterative methods lend themselves easily to the rntroductioij of 
additional correction fketers such as attenuation. With reference to FIGURE 1, the basic 
concept behind iterative methods, such as maximum likelihood expectation maximization 
(ML-EM) reconstruction, Is to find the object (La, image or emission gup) fbr wifich a 
mathematical projection produces results comparable to the set of measured emission 
projections. The ML-EM reconstruction algorithm is one such engine that allows the image 
to be found in an efficient manner. From an initial guess or assumption, a set of projec tions 
is artificially created by a projector employing a projection model. These projections are 
then compared to the "real" or measured set When certain conditions are mat (Le. a whan 
the projections of the emission map are sufficiently close to the measured emission 
projections), the iterative process stops and the cuxrent image is the best possible 
representation of the object, otherwise, the initial guess is updated and a new set of 
projections produced. Clearly, the projection model employed is an important part of the 
projection operation. However, an accurate projection model can be complex and detailed. 
[0004] 

With reference to FIGURE 2, one way to improve the "realism" of the projection model or 
the update operations is to include a priori information, Le,, what is already known about 
the object (e.g., hs contour, texture, and so forth). In general, this is a powerful tool that 
offers dramatic improvement of the reconstruction. One caveat is, however, that this 
info rm a tio n needs to be real In other words, if a priori information » true or accurate, it 
helps the reconstruction. On the other hand, inaccuracy tends to inappropriately bias the 
reconstruction or otherwise introduce unwanted artifacts. Moreover, in practice, it is 
difficult to find non-trivial characteristics of the object to be imaged that are always 
accurate. 
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[0005] 

With reference to FIGURE 3, the attenuation map is certainly one element that can hdp the 
reconstruction. Information pertaining to attenuation characteristics can be generated 
artificially for simple situations (e,g, uniform attenuation, symmetrical attenuation, etc). 
In general, however, inaccuracies can be introduced by suggesting attenuation features that 
are sot true and/or not object specific. 
[0006] 

With reference to FIGURE 4, it is significantly more advantageous when the at t ecnefl o n 
map is constructed for each object and used directly in the reconstruction- Often, the 
attenuation map is derived from measured transmission projection data. This is the basis , 
of most of the modern nan-uniform attenuation con ection devices. However, the problem 
with this approach if statistics. Typically, in order to define a usefal and usable attenuation 
map, many counts are needed. In practice, however, the definition of the at t e nu a t i o n map 
is count-limited, and the inherent noise associated with it is transported into the 
reconstruction of the emission map. Another problem is that often the transmission 
projection data used to derive the attenuation map is truncated thereby resulting in a 
degraded attenuation map. In other words, a poor attenuation map can actually degrade the 
image it aims at improving because of a lack of counts and/or a truncated transmission data 
set 

[0007] 

These limitations in the acquisition of the measured attenuation maps have provided the 
motivation for investigation into different techniques. It is generally known that the 
emission data contains information from which theoretically the attenuation map could be 
reconstructed. Such reconstruction methods are referred to as sourcetesa. In Act, each 
point in the object can be considered as a transmission source, and the observed intensity 
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at any given point on the detector can be compared v&h the expected intensity without 
attenuation. Unfortrmaarfy, due lWtsffans„ »it?h m th* IttmtH r m^ ^r of cou n ts in th» 
emission map, the relationship between the emission and attenuation maps can only 
partially be established. Moreover, then is an intrinsic mabiEty of the process to 
differentiate between a low activity, low attenuation condition and a high activity, high 
attenuation condition 

[0008] 

Yet another drawback of previous techniques is that they tend to employ only a linear 

projection model That is to say, they merely account for those events that fie along a 

s mg n l t r path or line from the detection bin to the activity. Such projection models are not 

adapted to account for scatter ox collimator resolution which may a detector 

bin to detect photons from off-line sources. 
[0009] 

[Problems that the Invention is to solve] 

The present invention contemplates a new technique for providing attenuation correction 

in emission computed tomography. 
[0010] 

[Means for solving the Problems] 

In accordance with one aspect of the present invention, a method of ML-EM image 
reconstruction is provided for use in connection 'with a diagnostic imaging apparatus that 
generates projection data. The method includes collecting projection data including 
measured emission projection data. Initial emission and attenuation maps are assumed and 
iteratxvdy updated. With each iteration, the emission map is recalculated by taking a 
previous emission map and adjusting it based upon: Q) die measured amission projection 
data; (U) * reprqjection of the previous emission map which is carried out with a multi- 
dimensional projection model; and, 02) a reprojection of the attenuation map. Also, with 
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each iteration, the attenuation map i* recalculated by taking a previous attenuation map and 
adjusting it based upon: Q) the measured emus on projection data; tod, (n) the reprojection 
of the previous emission map which is carded out with the muM-dimensional projection 
model 
[0011] 

In accordance with another aspect of the present invention, a medical diagnostic imaging 
apparatus for reconstructing image representations of a subject being examined therewith 
is provided. It includes an emission memory which stores measured emission projection 
data from the subject An image memory stores emission maps, and an attenuation map 
memory stores attenuation maps. A first projector generates emission map projections via 
forward projection of die emission maps from the image memory. Likewise, a second 
projector generates attenuation map projections via forward projection of the attenuation 
maps from the attenuation map memory. The first and second projectors use first and 
second iuute-c5meniional projection models, respectively. A first data processor samples 
the emission memory, and the first projector, and In accordance therewith updates the 
attenuation maps stored in the attenuation map memory. A second data processor samples 
the emission memory, the first projector, and the second projector, and in accordance 
therewith updates the emission maps stored in the image memory. Finally, a human- 
viewable display renders reconstructed image representations of the subject from the 

emission maps. 
[0012] 

[Embodiments] 

Ways of carrying out the invention win nw be described in detail, by way of example, with 
reference to the accompanying drawings. 
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With reference to FIGURE 5, the basic flow of a souroeless ML-EM reeoiistruetion Is 
accordance with the present invention is illustrated. A medical diagnostic imaging 
apparatus 10 (described in greater detail below) generates emission projection data. Inthia 
embodiment, the transmission map is reconstructed from emission data only. Lack of a 
transmission scan reduces the patient's radiation exposure. Moreover, in a preferred 
embodiment a multi-dimensional projection model provides a more accurate reprojection 
of the emission map. With reference to FIGURE 6, the basic flow of a source-assisted 
reconstruction is illustrated which is a flirther adaptation of the aourceless reconst ru cti on , 
In this embodiment, the medical diagnostic imaging apparatus 10 also generates 
transmission projections. Source-assisted attenuation correction combines the benefit of 
the real transmission measurements to the stability of the analytical sourceless approach. 
In this approach, the measured transmission projections, however statistically poor and/or 
truncated, are used to stabilize the analytics! approach and force the attenuation map and 
the emission map to be independent The transmission information and emission counts 
are combined to form an image of the attenuation map. The result is improved quality in 
the attenuation map. because the emission data is used to "assist" in its reconstruction. 
This, in turn, improves the reconstruction of the emission map in terms of quantitative 

accuracy. 
[00131 

The medical diagnostic imaging apparatus 10 u geoeralty any iwclear medicme aeanaer that 
produces scintigraphic images, for example, SPECT scanners or PET scanners. An 
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appropriate device is one that detects and records die spatial, temporal, tad/or otter 
characteristics of emitted photons. 
[0014] 

More specifically, with reference to FIGURE 7, in an exemplaiy embodiment, the 
^agnostic rudear Imaging apparatus or tcarmrr 10 is a gamma camera or SPECT scanner 
including a subject support 12, such as a table or couch, which supports and positions a 
subject being examined and/or imaged, such as a ph»ntnm or patient, within a subject* 
receiving aperture 13. The subject is injected with one or more radiopharmaceutical or 
radioisotopes such that emission radiation is emitted therefrom. A fint gantry 14 holds ■ 
rotating gantry 16 mounted thereto. The rotating gantry 16 defines the subject-receiving 
aperture IS. One or more detector heads 20 sre adjustably mounted to the rotating gantry 
16 with varying degrees of freedom of movtmeot Being mounted to the rotating gantry 16, 
the detector heads 20 rotate about the subject-receiving aperture 18 (and the subject when 
located therein) along with the rotation of the rotating gantry 16. In operation, the detector 
heads 20 are rotated or indexed around the subject to monitor radiation from a plurality of 
directions to obtain a plurality of different angular views. 
C0016] % 

Each of the detector heads 20 has a ratEaiion-receiviog fkee facing the aubjecfc*eceMng 
aperture II that includes a adntflladon crystal, such as a large doped sodium iodide crystal, 
that emits a flash of light or photons in response to incident radiation. An array of 
photomultiplier tubes receives the light and converts it into electrical signals. A revolver 
circuit resolves the x, y-coordinates of each flash of Kght and the energy of the incident 
radiation. That is to say T radiation ***** the scintillation «yw^T cmntog thff ftrfntffl^on 
crystal to scintillate, Le., emit light photons in response to the radiation. The photons are 
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cfiractad toward the photnrmiltrpficr tube*. Relative outputs of the photonaJtlpIler tubci art 
pnxwedaixlcoirtcttdtogtnanitcinou^ signal indicative of (1) a posmoii coordin ate 
on the detector head art which each radiation event is received, end (2) an energy of each 
event The energy a used to differentiate between various types of radiation such as 
multiple emission radiation sources, stray and secondary emission radiation, transmission 
radiation, and to eliminate noije. An image representation is defined by the radiation data 
received at each coordinate. The radiation data is then reconstructed into an image 
repres entat ion of the region of interest 
[0016] 

Optionally, the detector heads 20 include mechanical collimators 24 removably mounted 
on the radiation receiving feces of the detector heads 20. The collimators 24 preferably 
include an array or grid of lead or otherwise radiation-absorbent vanes which restrict the 
detector heads 20 from receiving radiation not travelling along selected rays in accordance 
with the data type being collected (Le., parallel beans, fan beam, and/or cone beam). 
Alternately, for example, when operating in a coincidence mode, the collimators 24 may 
be omitted. 
[0017] 

On e or more radiation sources SO ate mounted across the subject-receiving ap erture 1 8 flora 
the detector heads 20. Optionally, they are mounted between the detector heads 20 or to the 
radiation receiving faces of opposing detector heads 20 such mat transmission nufiation 
from the radiation sources 30 is directed toward and received by corresponding detector 
heads 20 on an opposite side of the subject-receiving aperture 18. In a preferred 
rmhoriimrnt, the collimators 24 employed on the detector heads 20, in effect, collimate the 
transnriidon radiation. Hiat is to say, tbe collimators 24 restrict the detector heads 20 from 
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receiving those portions of transmission radiation not travelling along raye normal (for 
parallel beam configurations) to toe radiation receiving faces of the detector heads 20. 
Alternately, other collimation geometries are employed aod/or the coffimation may take 
place at the source. 
[0018] 

In a preferred embodiment, the radiation sources 30 are moving point sources each 
traversing the length of the respective detector heads 20 to which they correspond 
Alternately, line sources are employed such as, e g, thin steel tubes filled with 
radionuclides and scaled at their ends. In other embodiments, the radiation source 30 are 
bar sources, point sources, flat rectangular sources, disk sources, flood sources, a tube or 
vessel filled with radionuclides, or active radiation generators such as x-ray tubes. 
Alternately, one or more point sources of transmission radiation may be utilized. 
[0019] 

With reference again to FIGURE 7. the running of an Imaging operation includes an 
iterative ML-EM reconstruction technique wherein emission projection data is 
reconstructed via an image processor 100 into an image representation of the distribution 
of radioactive material in the patient Of course, the reconstruction technique implemented 
varies slightly in accordance with the types of radiation collected, the types of collimators 
u»ed fan, cone, parallel beam, and/or other modes), and the diagnostic imaging 
apparatus employed. In any case, enuision radiation from the patient is received, for 
•wnple, by the detector beads 20 of the SPBCT scanner, and measured emission projection 
data is generated. The measured emission projection data normally contains inaccuracies 
caused by varying absorption or attenuation characteristics of the patient's anatomy. 
Optionally, a transmission scan is also performed such that transmission radiation from one 
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or more of the tr ansmi s sion nrfisti on sources 30 is also received, for example by the 
detector heads 20 of the SPECT scanner, and measured transmission projection data is 
generated. However, the measured transmission projection data is often noisy and/or 
truncated Where a transmission scan is performed, a sorter 200 sorts the measured 
e miiMinri tnd transmission Flection data on the basis of thair relative energies. The data 
is stored in a measured projection view memory 110, more specifically, in a corresponding 
emission memory 112 sod transmission memory 114. 
[0020] 

In those instances where transmission projection data is collected, a pre~acan blank run it 
performed to collect a baseline or reference with which to evaluate and/or compare the 
measured transmission projection data against. In the blank run, the s cann er is used to 
p erfbrm a transmission scan without the subject present therein such that a blank set of 
transmission projection data is collected. In a preferred embodiment, this data is stored in 
a reference memory 120. 
[0021] 

In preparation for the first iteration of the reconstruction process, sn attenuation map 
memory 130 and an image memory 140 are initialized by loading them with assumed or 
first estimate! of the a ttenuation map and emission map respectively, The first estimate fhr 
the attenuation map is optionally characterized by a uniform attenuation value inside a 
predetermined contour which contains the subject and zero outside the contour. Likewise, 
the first estimate for the emission map is optionally characterized by a uniform value inside 
the contour and zero outside. Alternately, the availability of additional a priori information 
allows for more accurate first estimates. 
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[0022] 

With each iteration of the reconstruction process, a first projector 150 reprojects and/or 
forward projects the estimated emission map stored in the image memory 140 to obtain 
projection views thereof Likewise, a second projector 152 reprojecti and/or forward 
prqjects the atfrrasfion map stoned in the attenuation msp memory 13 0 to obtain pro} ection 
views thereof: The projections are performed using muM-dimenaionai projection models. 
That is to say. the mulri-dimensional projection models account for radiation from multiple 
directions, as opposed to merely linear models that only account for radiation travelling 
along a singular direction. In a preferred embodiment, the projection models a cc o un t for 
any one or more of the following: scatter, the geometric point response (resolution) 
associated with the imaging apparatus, randomness corrections which account for random 
simultaneous photon detection not resulting from the same annihilation event in PBT 
applications, and/or cross-talk due to changing photon energies. In at least one 
embexfimant, die projection models used by the first and second projectors 150 and 1S2 are 
not the same. Alternately, for some applications, identic^ pixgection mcdeb are engaloyc^ 

A firat data processor 16 0 updates the atfcmi ation maps in tho attenuation map memory 130 
with each iteration, and a second data processor 170 updates the emission maps in the 
image memory 140 with each iteration. 
[0023] 

The firat data processor 160 recalculates the attenuation map by taking a previously stored 
attenuation map from the attenuation map memory 130 and adjusting It baaed upon a 
number of Actors. In a preferred embodiment, employing so called sourcelass 
reconstruction, the fectoxi include: (0 the measured emission projection data from the 
emission memory 112; and, (U) * reprojection of the previous emission map stored lo the 
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image memo ry 140 which U rcprqiectfid via the projector 150. More specifically, the first 
data processor 160 samples the emission memory 112 and the output projection views from 
the projector 150 to complete each update of the attenuation map which if then reloaded 
into the attenuation map memory 130. For source! ess reconstruction, measured 
transmission projection data is not utilized. As such, it is preferred that transmitting 
projection data not be collected to thereby minimize radiation exposure and reduce 
mechanical complexity. In an alternate embodiment, employing so called source-assisted 
reconstruction, in addition to the above, measured transmission projection data (optionally 
truncated) is also incorporated with each update of the attenuation map. In this 
embodiment, the measured transmission projections are also sampled from the transmission 
memory 114 by the first data processor 160 Car each update or recalculation. In addition, 
baseline or r e fer e n ce transmission projections from the reference memory 120 are also 
taken into account. A detailed trrount of the first data processor's operation is given below 
in mote mathematical terms. 
[0024] 

In both sourcdess and source-assisted reconstruction, the second data processor 170 
updates the emission map based on the same ftctors. These actors include: (x) the 
measured emission projection data; (H) a reprojection of the previous emission map; and, 
ffi) areprrjection of the attenuation map. More specifically, with each iteration, the seoond 
data processor 170 samples the emission memory 112, the first projector ISO, and the 
second projector 152. Based upon the data collected, the second data processor 170 
recalculates the emission map and reloads the updated emission map back into the image 
memory 140. Again, a detailed account of file second data processor's operation is gives 
below m more mathematical terms. 
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[0025] 

Upon completion of the finil iteration, m video processor 300 extract! or otherwise receives 
data from the image memory 140 wHch represents a reconstruction of the radionuclide 
distribution within the region of interest The video processor 300 then formats h into 
image representations of the subject for viewing on a human-viewable display, such as a 
video monito r 310 or other appropriate rendering device. 
[0026] 

In more mathematical terms, the operation of the image processor 100 is as follows. I1m 
iterative ML»EMrecc^ It has 

a robustness in reconstruction that is difEcukibr ether rccon^ 

based on the solution of linear systems. An important feature of the ML-£M al gorithm is 
that given certain mild restrictions appropriate for nuclear medicine emission imaging; it 
is applicable to any linear system of equations, whether h be over- or under-determined, 
consistent or inconsistent. In any event, after a cumber of iteration, the ML-EM algorithm 
converges to a solution that is the most likely or the most probable, given the original data. 
[0027] 

One of the appealing qualities of this iterative algorithm is that a general probability model 
of detection is used. The underlying radionuclide distribution is continuous, but thie signal 
or projection is detected at a discrete number of bins, or sampling locations. Id what 
follows, Fj represents the total detected photons in the i* detector location (La, the 
measured emission projection data). There sreM detector locations. Alio, the radionuclide 
distribution is represented as a discrete collection of basis functions. The basis function for 
the emission distribution is the voxel F k represent the average of the continuous 
ladionucfide distribution over the k* image voxel There are K image voxels. The transfer 
mat ri x coefficient (Le., the projection model) is the probability that a photon emitted in 
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voxel k per unit time is detected in detector location i. As each measurement can be noisy, 
the noise associated with each measurement is denoted by i*. Generally however, enough 
U known about the imaging system to be able to accurately model the statistical properties 
of the noise. Hence, the total photon count at the i* detector location is the sum of all 
voxels contributing to that, detector icea£oi3. This is mathematically expressed as: 

X 

where At, is the imaging time for thei* projection. 
[0028] 

This linear system of equations is then die basic imaging model Denoting fj** (Le., the 
k* voxel of the emission map) as the estimate of the activity F k in voxel k on the h* 
taxation of the MT^^reeonsto 
by: 

With this mediod scatter, resolution, snd attenuation can be incorporated into the model m^ 
However, the inclusion of attenuation and scatter In the transfer matrix employs a priori 
knowledge of the attenuation map. 
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[0029] 

Sonrcelei* Reconstruction for PET Embodiment 

The lack of <n attenuation map far emission imaging a this ease leads to a new formulation 
that now includes the unknown attenuation map. The forward problem is now stated as: 



where u, is the fineer attenuation coefficient of the j* pixel of the attenuation map, and c, 
ia the effective attenuation length of Uj. The let of pixels N M is given by. N„ » {pixel j | 
m, * 0, for source-detector pair Geometrically, the set N M is the set of pixels 

intersecting and covering the line of response i that contains the source pixel k. The 
estimate (f k ,u^of the unknowns (F k . Pj)" given by: 



- u*u _ f (*i » — 



P i m ki 

(4) , 



fhr the emission map which is updated via the second data processor 170; and 



for the atte nu ati on map which is updated via the first data processor 160. That is to say, in 
this embodiment (Le.. sourcetess reconstruction with PET data), the second data processor 
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170 and the first data processor 160 carry out the former tad later computations from 
equations (4) and (5), respectively. 

[0030] 

Source! ess Reconstruction for SPECT Embodiment 

For SPECT, m similar set of algorithms is developed. The forward problem in this case Is 
again: 



However, for SPBCT, the set N M is defined differently as illustrated in FIGURE 8. The act 



projection bin I This leads to a variation of the algorithm. In this embodiment, the 
emission update is given by; 




• • • • 



<«). 



N„ is tto act of fbcOs that cow the itteauttian lias integral fiom the •ourco pixel k to the 




and, the transmission image update is given by: 
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As before, in this embodiment (Le^ aourceless reconstruction with SPECT data, the second 
data processor 170 and the first data processor 16*0 cany out the former and later 
computations from equations (7) and (8), respectively. 
[0031] 

Source-Assisted Reconstruction for PET Embodiment 

Source-assisted reconstruction utilizes the measured transmission data, represented by T v 
and blank scan, represented by E^ to stabilize the sourceiess reconstruction. It 
accomplishes this by incorporating the transmission scan as a priori information in the 
likelihood probability Amotion^ The result is a different update fbr the transmission map 
only. For PET, the transmtwioo update is now given by; 

(9} , 

inwhich Atjis the scan time per projection for the emission data, and A* is a normalization 
constant that reflects the ratio of the emission scan-time to the transmission scan time. The 
coefficients £ yjt measure the contribution of the j* voxel of the attenuation map to than* 
transmission measurement and are not always the same as the coefficients defined above. 
They depend on the particular embodiment of the scanner employed. 
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10032] 

Source- Asiisied Reconstruction for SPECT Embodiment 

The up date rulo for the emission map in scurce-assisted SPBCT is the same as the up date 
rule for sourceless SPECT. Only die transmission map update rule is changed to the 
following: 

Optionally, the operationa of the first and second data processor* 1*0 and 170, as well as 

the other components, are implemented as hardware, software, or combinations of hardware 

and software configurations. 
[0033] 

[Advantages of the Invention] 

One advantage of the embodiments described above of source-assisted attenuation * 
correction for emission computed tomography, is improved attenuation maps for ML-EM 
reconstruction. Another advantage is reduced patient exposure to radiation via sourceless 
ML-EM reconstruction. Yet another advantage is rnulti-dxmensional projection modelling 
which more accurately approximates actual conditions. Another advantage is that there is 
less likelihood of noise being introduced into the final reconstruction via the attenuation 
map. A further advantage is that acc ur at e reconstruction is possible oven with truncated 
transmission data. 
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[0034] 

[Brief Description of the Drawings] 

[FIG. 1] is t flow diagram Illustrating the operation of a typical ML-EM reconstruction 
method in accordance with prior art techniques; 

[FIG. 2] is a flow diagram illustrating the operation of a typical ML-EM reconstruction 
method employing a priori information in accordance with prior art techniques; 

[FIG. 3] is a flow diagram illustrating the operation of a typical ML-EM reconstruct^ 
method employing an artificially generated attenuation map in accordance with prior art 
techniques; 

[FIG. 4] is a flow diagram illustrating the operation of a typical ML-EM reconstruction 
method employing an a ttenu a ti on map derived from measured t ransmissi on projection data 
in accordance with prior art techniques; 

[FIG. 6] is a flow diagram illustrating sourceless ML-EM reconstruction in accordance 
with aspects of the present invention; 

[FIG. 6]) is a flow diagram illustrating source-assisted ML-EM reconstruction in 
accordance with aspects of the present invention; 

[FIG. 7] is a d i a gr a mm a t i c illustration of a diagnostic nuclear im aging apparatus 
emric^g an image prccesra 
aspects of the present invention; and 

[FIG. 8] is a graphical flhutratian showing the set of pixels employed in the 
reconstruction of SPBCT data in accordance with aspects of the present invention. 



(35) 




4*B3 2 000-180550 




(37) 



$3rgB2 000-180550 




(40) 



$568 2000-180550 




(42) 



$$Sa 2 000-180550 




FIG. 8 
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[Document Namp] ABSTRACT 

[Abstract] 

[Problems] 

There is provided a new technique for providing attenuation correction in 
emission computed tomography. 
[Solving Means] 

A method of ML-EM image reconstruction is provided for use in connection with a 
diagnostic imaging apparatus (10) that generates projection data, The method includes 
collecting projection data, including me a sur ed emission projection data, An initial emission 
map and attenuation map are assumed. The emission map and the attenuation map are 
iterativdy updated. With each iteration, the emission map is recalculated by taking a 
previous emission map and adjusting it based upon: 0 the measured emission projection 
data; 00 a "projection of the previous emission map which is carried out with a multi- 
dimensional projection model; and, (Hi) a projection of the attenuation map. As well, 
with each iteration, the attenustion map is xecslnifirfM by taking a previous attenuation map 
and adjusting it based upon: ® the measured emission projection data; and. (H) a 
repiojection of the previous emission map which is carried out with the rmilti-dimcnsional 
projection modeL In a preferred embodiment, with source-assisted reconstruction, the 
recalculation of the attenuation map is additionally based upon: (Hi) measured transmission 
projection data; and <jv) a reference or blank data set of measured transmission projection 
data taken without the subject present in the imaging apparatus (10). 
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